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ABSTRACT 

Proper motions for several Milky Way dwarf spheroidal (dSph) galaxies have been determined 
using both ground and space-based imaging. These measurements require long baselines and repeat 
observations and typical errors are of order ten milli-arcseconds per century. In this paper, we utilize 
the effect of "perspective rotation" to show that systematic proper motion of some dSphs can be 
determined to a similar precision using only stellar line-of-sight velocities. We show that including 
the effects of small intrinsic rotation in dSphs increases the proper motion errors by about a factor 
of two. We provide error projections for future data sets, and show that proposed thirty meter class 
telescopes will measure the proper motion of a few dSphs with milli-arcsecond per century precision. 
Subject headings: Cosmology: dark matter, theory-galaxies: kinematics and dynamics-Astrometry 
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1. INTRODUCTION 

The three-dimensional motions of stars and galaxies 
provide valuable information on the local Universe, rang- 
ing from planetary companions of nearby stars to the 
orbital properties of near by galaxies and Galactic satel- 
lites (jUnwin et al.l 120081 ). However, because measure- 
ments of proper motion must be done with respect to 
cosmic standards of rest, such as background galaxies 
and quasi-stellar objects, they depend sensitively on the 
number and nature of background objects in the target 
field. To obtain measurements of required precision, long 
baselines and repeat measurements are necessary, rang- 
ing from several years for space-based telescopes to tens 
of years for ground-based telescopes. 

Nonetheless, despite their low luminosities, large 
distances, and small angular separations, about a 
dozen nearby galaxies n ow have measured proper mo- 
tions (jPiatek et al.ll2007l ). Among the nearest of these 
galaxies are the dwarf spheroidals, which have observed 
luminosities that vary anywhere from a thousand to a 
million times the luminosity of the Sun. These galax- 
ies are supported primarily by their velocity disper- 
sion and have high mass-to-light ratios ^ 100 Mq/ 
Lp, fMateo"1998': 'Lokas et al."2G05':'Gilmore et al .ll2007l: 
Strigari et al. 2007; Simon & Geha 2007; Walk er et all 
20071) . Proper mo tions have been measured for s everal 
Milky Way dSphs (iPiatek et al.l[2(M l2?M l200l . with 
errors ~ 10 milli-arcseconds per century, corresponding 
to transverse velocity errors ^ 100 km s~^ for typical 
dSph distances. In contrast, their velocities in the direc- 
tion of the observer are now known to ^ 1 km s~^. 

In this paper, we introduce an alternative technique 
for determining the proper motions of Milky Way dSphs. 
Specifically, we utilize present samples of stellar line- 
of-sight velocities together with the effect known as 
"perspective rotation," in which the tangential mo- 
tion of the galaxy contributes to the measured line-of- 
sight velocity at large angular separations from the cen- 
ter of the galaxy. Perspective ro tation has been de - 
tected in Galactic globular clusters (Me rritt et al.lfl997f ) , 
and has been used to measure the distance to the 
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Large Magell anic Cloud iG ould' '2000") and the mass 
of M31 ( van d er Marel fc Gu hathakurta 2007). Here we 
show that, using perspective rotation, proper motions of 
several dSphs can be determined to a precision rivaling 
the best existing ground and space-based measurements. 
Additionally, the first proper motion measurements will 
be possible for several dSphs that are difficult to access 
via traditional methods. 

2. PERSPECTIVE ROTATION 

Perspective rotation is simple to understand if we con- 
sider the dSphs as extended objects. Because of their 
close proximity and spatial extent, the line-of-sight ve- 
locities of the stars vary as a function of projected ra- 
dial separation, R, from the center of the galaxy. As 
R increases, the line-of-sight velocities receive increas- 
ingly larger contributions from the tangential motion of 
the object in space. The net result of this radially vary- 
ing li ne-of-sight ve locity is known as perspective rotation 
(jFeast et al.lll96i[ ). An object that is not rotating intrin- 
sically acquires a velocity gradient that is proportional 
to the projected distance from the center of the galaxy. 

To describe the effect of perspective rotation, we define 
a cartesian coordinate system in which the z-axis points 
in the direction of the observer from the center of the 
galaxy, the a;-axis points in the direction of decreasing 
right ascension, and the y-axis points in the direction of 
increasing declination. The angle (j) is measured counter- 
clockwise from the positive x-axis, and p is the angular 
separation from the center of the galaxy. The linc-of- 
sight velocity is then 

vios = Vx sin p cos (p + Vy sin p sin (j) — cos p. (1) 

For all of the dSphs that we consider, it is appropri- 
ate to use the small angle approximation, sinp ~ R/D, 
where R = ya^^+j/^, and D is the distance to the dSph. 
Then using smcf) = y/R, equation [1] can be written as 
vios = Vxx/D -\- Vyy/D — Vz- In the limit that sinp ^ 1, 
the line-of-sight velocity is constant across the dSphs and 
we recover in this limit that vios = —Vz- It is evident from 
equation [T] that the transverse velocities and Vy have 
the maximal contributions for galaxies that are the ap- 
propriate combination of the most nearby and the most 
spatially extended. 
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3. LIKELIHOOD FUNCTION AND ERROR PROJECTIONS 

We define the Hkehhood function for an observed set 
of line-of-sight velocities as 

N ^ r 1 f ^2^ 

2 



m = n 



1 



=1 



: exp 



(2) 



Here 6 is the set of parameters that describe the model 
of the galaxy. The product is over the total number 
of stars, with measured velocities, and i is an index 
that represents a star that is located at a fixed projected 
position. The total velocity dispersion, a, is the sum of 
the intrinsic dispersion, aios, and the dispersion from the 
measurement, am- 

We determine aios from standard dynamical equilib- 
rium analysis, assuming that the potentials of the dSphs 
are spherically-symmetric. The Jeans equation for the 
radial velocity dispersion cr,. is 



,GAf(r) 

r — 1 = -P*[r) 

ar r 



2P{r)p4r)al (3) 



Here M{r) is the halo mass, (3{r) — 1 — (Jg/a'^ is the 
stellar velocity anisotropy, and p^, (r) is the three- dimen- 
sional density for the stars, which is determined from 
the measured surface density of stars, I{p). For I{p), 
we use King profiles, which are good fits to the sur- 
face densities of the dSphs that we examine (|Kindll962l : 
llrwin fc Hatzidimitrioul[l995f ). King profiles are fit by a 
"core" radius, Tc, and a limiting radius. Him. For a King 
profile the limiting radius is the tidal radius of the stars. 

The radial velocity dispersion can be determined by 
solving for ar using equation [3] and imposing the bound- 
ary condition (j^ — > as r ^ oo. We can then use 
equation [3] and the definition of f3 to convert the veloc- 
ity dispersions implied from equation [1] into observable 
quantities by integrating along the line-of-sight through 
the dSph. Performing the integration gives 

2 f°° „, . oi Pi,{r)a^rdr 



Hp) 
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= cos^[6'(r)] sin^ p + 
sin^ [9{r)] cos^ p, with 



Here, rmin = I?tanp; i 
2 sin[0(r)] sin p cos[9{r)] cos p 

cos[6{r)] ~ {D — D cos^ p+ \/r'^ ~ sin^ p cos p)/r. For 
simplicity we assume a constant /3; we find that our re- 
sults below do not depend on whether /3 is a fixed con- 
stant or is allowed to be a more complicated function of 
radius. Accounting for the definition of the line-of-sight 
velocity in equation [1] equation 2] differs from the stan- 
dard def inition of the projected ve locity dispersion at a 
fixed R (jBinnev fc TremaindfTOSTi) . 

We take the dSphs to be dark matter dominated, 
with halos described by the density profile p(r) — 
Poir/royil + ^Vr'[;)("~'='/^ We take the slopes, a, b, 
and c, the scale density po, and ro to be unknown pa- 
rameters and marginalize over them with uniform pri- 
ors. Our results are independent of the halo mass model, 
provided the projected velocity dispersion is fixed to 
match the observed, nearly flat dispersion profiles of the 
dSphs (I Walker et al.l[2007li . 

We are interested in using the likelihood function in 
equation [2] in concert with equations [3] and |4] to project 
the errors on the components of the transverse veloc- 
ity, Vx and Vy. The attainable errors depend on the 



covariance matrix for the model parameters 6, which 
we will approximate by the Fisher information matrix 
Fab = {d^ In C/d0adeb) (iKendall fc StuartI [19691) . Our 
model parameters are those that describe the dark mat- 
ter halo, the velocity anisotropy, and the spatial motion 
of the galaxy. The inverse of the Fisher matrix, F^^, 
provides an estimate of the covariance between the pa- 
rameters, and vi^Mi approximates the error on the pa- 
rameter 9a. The Cramer- Rao inequality guarantees that 

\/ Faa is the minimum possible variance on the a*^ pa- 
rameter for an unbiased estimator. Using F^^ in place 
of the true covariance matrix involves approximating the 
likelihood function of the parameters as Gaussian near 
its peak, so F~^ will be a good approximation to the 
errors on parameters that are well-constrained. 

We construct the Fisher matrix by differentiating the 
log of the likelihood function in cquation[21 and averaging 
over the data. Performing the appropriate averaging, and 
using the above definition of the line-of-sight velocity in 
equation [1] our final expression for the Fisher matrix is 
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In deriving equation [Sj we have assumed no correlations 
between the theory and the measurement dispersions. 

In the second term in equation [5l the derivatives are 
with respect to the theory dispersion alone, whereas both 
of the contributions to the variance sum in the denom- 
inator. For the classical, well-studied dSphs, the intrin- 
sic velocity dispersions are of order 5-15 km s~^, while 
the mean measurement uncertainty is less than 2 km 
s~^, so the dominant contribution to the variance comes 
from the theoretical distribution function (For many of 
the newly-discovered satellites, h owever, both contrib u- 
tions to the dispersion are similar ijSimon 
Here we are interested in the highest luminosity dSphs, 
so in all of our examples the error from the intrinsic dis- 
persion of the dSph is the dominant contribution. 

Equation [5] shows that, to determine the error on any 
of the 9 parameters, we need two pieces of information: 
1) the distribution of stars within the dSph that have 
measured velocities, and 2) the error on the velocity of 
each star. The projected errors are independent of the 
mean velocity of the stars. The errors on the param- 
eters describing the dark matter halo and the velocity 
anisotropy enter only through the second term in equa- 
tion [5] when we differentiate the velocity dispersions in 
equation m The errors on the three velocity components 
Vy, and Vz are independent of the parameters that 
describe the halo and velocity anisotropy. 

4. RESULTS 

To project errors on the proper motion components, 
we draw stars from uniform distributions for both the 
X and y spatial coordinates. We assign a measurement 
error of Um = 2 km s~^ for each star, which represents 
a conservative upper limit for the mean error in high 
luminosity dSphs. Our results arc not strongly sensitive 
to the shape of the distribution from which we draw stars, 
except in the case that the stars are strongly centrally- 
concentrated within the King core radius. 

We consider two fiducial models for dSphs. The first 
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model is a compact, nearby system designed to repre- 
sent Draco, which is located at a distance D = 80 kpc 
and is described by a King core and limiting radius of 
0.18 and 0.93 kpc, respectively. From ground based mea- 
surements, the proper motion of Drac o is (mchM^) = 
(60 ± 40, 110 ± 50) mas per century (jScholz fc IrwinI 
Il994f ). We note that Sculptor is at a distance similar 
to Draco, has a higher luminosity, and contains a similar 
number of line-of-sight velocities, so it may also present 
an interesting target. However, as we discuss below, 
the prospects for measuring its proper motion using our 
methods may be complicated b y the presence of intrinsic 
rotation (jBattagha et al.ll2008^ . 

As our second model, we consider a more distant and 
more extended dSph, designed to represent Fornax (D = 
138 kpc). Fornax is described by a King core radius and 
limiting radius of 0.40 and 2.7 kpc, respectively. The 
proper motion of Fornax from Hubble Space Telescope 
(HST) observati ons is (mc.M^) = (4 7.6±4.6, -36.0±4.1) 
mas century" -'^ (jPiatek et al.ll2007h . At present, Fornax 
has a total of ^ 2 000 measured line-of-sight velocities 
(| Walker et al.ll2007l ). and ^ 8600 red giant stars greater 
than 20th magnitude for which future observations will 
be possible. 

Examining equation [5l the projected errors on param- 
eters generally depend on the fiducial point in parameter 
space. We consider fiducial models fit to match the pro- 
jected velocity dispersions of Fornax and Draco. Wc find 
that po = lO"^ Mq kpc~^, tq = 2 kpc, a = 1, = 1, 
and c = 3 provides a good description of their disper- 
sion profiles (Strigari ct al. 2007). For Fornax we take 
the stellar distribution to be isotropic, (3 = 0, while for 
Draco we take a tangential anisotropy of (3 = —1. The 
dark matter halos have a mass of ~ 10^ Mq within their 
inner 300 pc. 

In Figure[T]we show the projected errors on the proper 
motion components as a function of the number of line- 
of-sight velocities. For each model, for numbers of stars 
^ 1000, the errors reduce to tens of milli-arcseconds 
per century, r ivaling the best measurements from HST 
i|Piatek et al l [20 5, 2006 . 2007i ). For the specific case of 
Draco, we proj ect that the curren t samples of line-of- 
sight velocities (jWalker et al.ll2007f ) will reduce the error 
on its proper motion by a factor of at least two relative 
to the best ground-based measurements. For each curve 
in Figure [1] we show only one component of the proper 
motion, in this case the component corresponding to the 
transverse velocity Vx, and note that the constraints on 
the component corresponding to Vy are similar by sym- 
metry. We also find Vz to be strongly constrained, <^ 1 
km s~^ for the samples of stars that we consider. 

It is important to note that the projected errors in Fig- 
ure [1] are valid only for the spatial distributions of stars 
and the measurement errors that we consider. The true 
constraints will depend on the detailed distributions of 
these quantities, and may in fact be reduced with more 
precise velocity measurements or varying stellar spatial 
distributions. For example, using the sample of Fornax 
red giants greater than 20th magnitude, we find sensitiv- 
ities reduce to ~ 3 mas per century. 

5. CONTRIBUTION FROM INTRINSIC ROTATION 

In the analysis above, we have assumed a negligible 
contribution from the intrinsic rotation in dSphs. This 
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Fig. 1. — The projected errors on the proper motion of dSphs as a 
function of the number of line-of-sight velocities, for two different 
models. The left panel uses the Draco fiducial model, while the 
right panel uses the Fornax fiducial model. In each panel, the solid 
(red) curve is for a non-rotating system, and the dashed (blue) 
curve includes a small rotational component. The lines indicate the 
approximate present errors on the proper motions. Measurement 
errors of 2 km s~^ are assumed for line-of-sight velocities. 

is a good model for the majority of these systems, which 
exhibit no detectab le rotational or streaming motions in 
the kinematic data (jWalker et al.ll2006l:lKocli et al.ll2007f) 
(See however the recent results of IBattaglia et al.l ([2001), 
which show that rotation may be present in Sculptor). 

Though the intrinsic rotation in dSphs is small, it is 
important to determine how even a small signal may de- 
grade the constraints on the proper motions we have de- 
termined above. Probably the simplest rotational model 
to consider is a sinusoidal variation in the rotatio nal ve- 
locity as a function of azimuthal angle (Drukier et al.l 
Il998l) . In this model, a term of the form ylsin(0i -|- 00 ) 
can then be added to equation [l] where A is the ampli- 
tude of the rotational motion, and is the projected axis 
about which the rotation occurs. We note that higher or- 
der multipolcs may exist in the velocity field as a result of 
more complicated rotation or tidal effects; however here 
we assume these higher order terms are sub-dominant to 
the leading dipole term. 

We add the parameters A and (po to the Fisher ma- 
trix, and marginalize over them with uniform priors. In 
Figure [TJ we show the projected proper motion errors 
including the effect of intrinsic rotation. Accounting for 
intrinsic rotation, we find that the errors on the proper 
motions may increase by about a factor of two for a fixed 
number of stars. Again, the details of the constraints de- 
pend on the exact distribution of the stars. Regarding 
the rotational parameters themselves, we find that is 
not well-constrained, but A is determined to a precision 
of ^ 1 km s""'^ with ^ 500 line-of-sight velocities. Sim- 
ilar to Vx and Vy, these errors on A are independent of 
its mean value. We note that for Draco and Fornax the 
data suggest A <1 km s~^. 

We reiterate that our parametrization of rotation is a 
toy model introduced to understand the effect of includ- 
ing rotation on the precision with which proper motion 
may be measured. Although we have determined the 
errors on A by solving for aios using the Jeans equa- 
tion, this procedure is not self-consistent, because in the 
presence of rotation the Jeans equation itself must be 
modified. In a rotating system, this would be impor- 
tant for determining the parameters describing the halo 
or the velocity anisotropy. However, since here we are 
interested only in determining the velocity components, 
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all that we demand is aios have enough freedom at each 
projected position to fit the kinematic data. As long 
as the streaming motion is small, as in the case of the 
dSphs we consider, this procedure should provide a good 
approximation to the errors on the model parameters Vx, 



and A. 



6. CONCLUSIONS 



We have shown that line-of-sight velocities from dwarf 
spheroidals can determine the transverse velocities of 
these galaxies to a precision of ^ 100 km s'-^. This mea- 
surement utilizes "perspective rotation," or the variation 
of the line-of-sight velocity across the galaxy resulting 
from its proper motion. The above sensitivity is simi- 
lar to the current measurements from HST and ground- 
based imaging for dSphs such as Fornax, Carina, and 
Sculptor. Using perspective rotation, the proper motion 
errors for several dSphs, including Draco, will be reduced 
by at least a factor of two, and proper motions for dSphs 
such as Sextans, Leo I, and Leo II will be determined 
for the first time. We find that, with ~ 10^ stars, the 
proper motion of Sextans can be determined to a sensi- 
tivity of 5 mas per century, while for Leo I and Leo II 
we project errors of 10 mas per century with 10"^ stars. 

Prospects are promising for improving on these mea- 
surements with future, larger samples of line-of-sight 
velocities. For example, we project that with ^ 10* 
line-of-sight velocities from Fornax, the errors are re- 
duced by ^ 40%. Proposed thirty meter class tele- 
scopes, such as the Thirty Meter Telescope (TMT) 



(lhttp://www. tmt.org[) and the G iant Magellan Telescope 
(GMT) ( http: / / www.gmto.org[ ) , may measure veloci- 
ties for lO"* stars in multiple dSphs, obtaining milli- 
arcsecond per century sensitivity. Nearby ultra-faint 
dwarfs, with extent ~ 100 pc and distances 50 kpc 
(Coma Berenices, Willman 1, Ursa Major II), may also 
be studied. In these objects we project errors ~ 100 mas 
per century with the optimistic scenario of 10^ stars. 

A determination of dSph proper motions will provide 
a full three-dimensional mapping of their orbits, with 
prospects of improving the measurement of the Milk y 
Way dark matter halo mass (|Little fc Tremain^ 119871 ). 
In addition, the proper motions will allow for a detailed 
comparison to the orbital propert ies of dark matter sub - 
halos in numerical simulations (|Diemand et al.l l2007f ). 
and provide information on the merger and accretion 
histories of these satellites. It will be possible to de- 
termine which of the satellites are bound to the Milky 
Way, a property that may be closely intertwined with 
the complex star f ormation histories of several satellites 
(|Besla et al.l 120071: iMateo et al.ll2008[ ). 
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